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Abstract 

A solution for the mathematical problem of functional calculus 
with Laplace-Beltrami operator on surfaces with axial symmetry is 
found. A quantitative analysis of the spectrum is presented. 

1 Introduction 

The physical situation which has initiated this research is that of a dielectric 
particle with electric charges on its surface, placed in electric field. Here, 
the diffusion equation of the charges is coupled with the Maxwell equations. 
There is an analytical solution of this system of equation |IJ which involves 
some functional calculus with operators, in particular with Laplace-Beltrami 
operator defined on the surface of that particle. We can imagine many other 
physical situations described by a complicated system of equations where 
the Laplace-Beltrami operator is implicated (e.g. that of the acoustic wave 
scattering on particles with membrane, etc.). As before, one can find a 
compact solution by using functional calculus. However, these solutions are 
not complete because, at this level, all is formal. We must have an effective 
procedure to calculate the expressions which involve operators. One can try 
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to compute the matrices of those operators in some orthonormal basis and 
to transform the problem in to algebraic one. The practical problem is that 
one can compute only a finite number of matrix elements and this can leads 
to serious problem when unbounded operators are implicated. If we choose 
an inappropriate basis, it is possible that the expressions, calculated with 
truncated matrices, to not converge at the correct result. 

In this paper we will find an orthonormal basis in the space of square 
integrable functions defined on a surfaces with axial symmetry such that 
the truncated matrices of Laplace-Beltrami operator converge in the norm 
resolvent sense. Then, according to || we can use these truncated matrices 
in functional calculus. 



2 The result 

Let M be a C°° closed 2-dimensional surface which in the spherical coor- 
dinates {r, 9, <f)} relative to a 3-orthogonal system of axes is defined by the 
equation r = r (9). We consider that all necessary conditions to have a C°° 
surface are fulfilled. Let this surface be equipped with the metric which is 
induced by the embedding in R 3 and let x G M be the point defined by 
9 = 0. Relative to this point, the normal coordinates {A, tp} are defined by 

I A (x (9, </>)) = d{x (9, <f>) ,x )=J dt^J r {tf + r' (t) 2 

I <P = <P 

which parameterize the entire surface, without the points 9 = 0,tt. We 
define R = X(x (tt)) /it and the new coordinates: = X/R, if}. In these 
coordinates, the metric form is 

9 ^ ,(f ^ = ( r(#(tf)) 2 sin(#(tf)) 2 ) • 



Proposition 1 The set of C°° functions: 
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is an orthonormal basis in L2 (M, /j, g ), where Y\ m represent the spherical har- 
monics and fx g is the measure induced on M by the metric g. 



Proof. The orthonormality: 

n 2tt 



Mm, y Vm >) =j dd J d v v^et^ • y lm (■&, if) y*, m , if) = 





7T 2"7T 

J d§ J dpRr(9) sin# 



r sin 9 R R 





TT 2lT 



J dd J dp nn#-yi m (&,<p)y;, rn ,(&,<p) = 5 w 8 ri 





The completeness: 

it 2ir 



(] l,m 



TT 2TT 

J d& J dp' Rr {9') sin {&) ^ 

l ' m 



Rsiwd / Rsiwd' 

■x 



r(9) sinfl V r (9') sinfl' 



WW) ffl , n _Rsintf 

— r R—fV'V^ym^e* 

TT 2lT I 

I do'/ V E ^ (*» ^ W> f'nl%^ V) = f V, <p) , 







because yftjgfifif (0, </>) is in L 2 (M, if / e L 2 (M, N ). ■ 

For a fixed m, let S m be the Hilbert subspace spanned by {yim}i>\ m \, 
which is invarianted by the Laplace Beltrami operator. In the following, we 
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will consider the restriction of this operator at a S m subspace, A^"^ = A |s m . 
Let P^ m \ k > \m\, be the projection on the subspace spanned by the vectors 
y\ m \m, •••,34m- Our main result is: 



Theorem 2 The sequence of operators 



,(m) 



p(m) oA (ra) Q p(m) 



k>\m\ 



converges in norm topology at the operator [A*™) — z\ 1 , for any z G C with 
Imz 7^ 0. 



Proof. We have successively: 



,(m) 



pM A (m) pM _ ^ 



A( m ) — ^ 



P 



(m) 



pM A (m) pM _ 2 



A( m ) — ^ 



_ ( j _ pM 



1 



A( m ) - ^ 



* pM A (m) _ p(m)\ 



A (m) _ 3 i 7 ^ 



(m) \ 1 
' c 7 A (m) _ 



pM A (m)p(" 



r i I pHA(m) A r _ pW\ f j _ p{m)\ 1 

1 + z k V ^ ) V ^ JaM-z 



Without loss of generality we can choose z = iou, uj e R, u> 7^ 0. Thus: 



,(m) 



fc p[ m) o A(™) o P< m) - z AM 



< 



' 1 



)M A (m) ^ 



/-P. 



(m) 



/-P 



>(™A L_ 

* 7 A (m) _ 
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Lemma 3 For I, k E N, I, k > \m\, with I ^ k: 

(y lm ,A ( ^y km ) = (y lm ,h-y km ), 

where h is at least a C° function on M. The operators P^ o A*™) o 



I-P, 



(m) 



are bounded and their norms satisfy: 



pjjn) Q A ( m ) (j D (m) 



< 



It follows that 



P, 



pW o A( m ) o P ( k m) - z ^ {m) - z 



< 



1 + 



CO 



I-P 



(m) 



1 



A( m ) - i 



iuj 



To evaluate the last norm, we use the following 



Lemma 4 Let s{$) be the quantity ,Jjg^. //Ag,..., \ { ™ ] are the 

ordered eigenvalues of A^ and v 1 ,™) , . . . , v^j.-.are the corresponding eigen- 
vectors, then for any m e Z anrf / > \m\, I > 1 and n > \m\ 

Rc 



\(y lm \v^)\< 



IHIoc + c V A 



(m) 



and for any m e Z, Z > |m|, ano? n > \m\, n > 1 

II^IL + lv 7 ^"^) 



|<^K m) >|< 



A 



(to) 



Moreover: 



+ < A (m) < c 2 [j[+l) 



i? 2 



i? 2 



w/iere c = [max {p 2 ^ , ||s : 



}] 



1/2 
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Now, let v G L 2 (M, fi g ), v — ^2 a n ■ v n m ^ ■ Then 

n>|m| 



E E a ' 

l>k+l n>\rn\ 



yim,V: 



(m) 

'n 



A„ — tOO 



E 

i>k+\ 



E a * 

n>\m\ 



yim,Vn 



(m) 



\ ( m ) 

An — IUJ 



< 



E E i°«i 2 - E 

l>k+l n>\m\ n>\m\ 



yirai Vn 



\ ( m ) 



< 



n 2 E E 

l>k+l n>\m\ 



cR 



\\ds\L + c V X 



*> 2 !tt e 



n>|m| 



IHIoo + C\A 



. (m) 



An — ^ 



For \m\ > 0, 



< 



2 M 



E 



c + 



||<fe|l 



V 7 ^ 



fc+1 t^, Ai m) 1 + 



< 



n>|m| An x 1 A<T )2 



k + 1 



y — 

, (m) 
An 



n>|m| 
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Finally 



(to) 



I-p, 



For \m\ = 0, 



J AM - 



iu> 



< 



(cRf 



\m\ 



c + 



cR\\ds\ 



y/\m\(\m\ + l) 



( 



fc + 1 



(cR) 2 \\ds\t 
u 2 



< 



c + 



\\ds\\ c 



n>l A ™ i 



< 



,(0)2 



thus: 



(cfl)2 JM_ l!i + (ea) . 



c + 



^ ^ ) A (o) - 



< 



v^TT) V (ci?) 



2 2 

or 



+ 



c + 



cR \\ds\ 



Having that Hds^ = ^ ||ff = ^ ll s 'lloo we can conclude: 



3 M * 1 

k P fc (m) o AM o P fc (m) - icj AM - iu 



< 



1 + 



(cR) 2 



R»° »oo i ^2 



1 + 



V2 



, for m = 



1 + 



i + 



, for |m| > 1 



Proof of Lemma 3. We have successively 
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7T 2w 

J d'd J dip R 2 sin^s 2 



7r 2n 



1 d fY^\ d (Y km 



R 2 d'd V sR J d'd V sR 



m s 



2„- 6 y 

1 lm 1 km 



R 2 sm$ 2 R R 



/ dp- 



sin i? 







dd lm dr km 



din s d 



1 

R 2 



it 2tt 

d'd / dtp sin ($) 



2 — 4 

m s 



i? 2 sin 2 # 



km) 



Y* Y 

2 a I lm I km 



91ns 

dd 



Y* Y, 

1 Im 1 km 



+ 



7T 2lT 

J d'd J dipsmi!) 



2 — 4 

m s 



d + d 



+ 





7T 2TT 

— J d'd J d(fi sin $ 



7T 2lT 



dins 



1 d 

sin i? cfr? 



sin-$ 



dins 



' 1 Im 1 km 



dip 



sin-$ 



s - 4 -l 



m 



+ 



dins 



1 d 



sin-$ 



dins 



R 2 sin 2 $ ' \ d& J smtidti 

o o 

It is easy to check that s (i?) is at least of C 2 class, so the function h : [0, 7r 

R, 



* Im-tkm- 



h ($) = m 2 



1 



i? 2 sin 2 

is at least of C° class. Finally 



dins 



1 d 
sin (9i9 



sin-$ 



dins 



2w 



(y lm ,M m *>y km ) = -^ J dd J dct>sind-h{d)-Y^Y km 







2tt 



^ J dip sin (t?) s 2 • h (■&) • = (y lm , h ■ y kr , 



For the second part 

pM o A (m) o _ pWj 



I/, 
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'Pt\h\l-P^)u 



< 



r>( m ) 
P> >V 



h- I 



u 



< 



— II lloo 



\V • \\U\ 



Proof of Lemma 4- 



Proposition 5 The application g : (0, ir) x [0, 2n] — > M (2 x 2) 



9 <P) 



1 

i? 2 sin(tf) ; 



defines a metric on M. Moreover, y^lff = s2 - 



If we consider the spaces of the squared integrable functions with the mea- 
sures induced by the two metrics, L 2 (M, fx g ) and L 2 (M, fig), and the spaces 
of one-differential forms with the standard scalar products, (M, fi g ) and 
(M,fig), then: 



Proposition 6 The spaces L 2 (M, fi g ) and L 2 (M, fig) coincide, A^ (M,/i s ) 
and A^ (M, fig) coincide too. 



Proof. For / G L 2 (M, fx g ) we have: 

7T 2-7T 

^ / dyVdet^)imy?)| 2 < 



Vdet^(i?) 



Vdet^(i?) 



o o 

thus: f E L 2 (M, /i^) .Analogous, for f E L 2 (M, /^) results 

7T 2-7T 



=y y d<py/d^w\f(&,<p)f 



< 





thus f eL 2 (M,fi g ). 



VdrtgW) 



< oo, 
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Let uj G (M, fig), uj = uj$d<d + u^dtp. Will follow 

7T 27V 7T 2"7T 

u, to) = J d$ J d<p^/detg(<d) 





M + 



,2 

2 . I^VI 



det g 



< 



' det g / det g 



|2 



< max < 4 / , 4 / - — - > • \\u L < oo, 

5 J det 5' Vdet^J " " 9 ~ ' 

thus u; G (M, /^) . The same steps can be followed to show that uj G 
(M,fjig) uj G (M, Denoting 



c = \J max {ll s2 lloo - II s 2 lloo}- 



we have on L 2 (M, fx g ) = L 2 (M, /^) 

1 

c 

and, on A« (M, fi § ) = A™ (M, n g )\ 

1 

c 



9 — H H9 — c II Us 



9 — II H9 — c II Uj? 



Now, we have successively 

'Y, 



Jrn 1 ( m ) 
Si? 1 n 



R 





' 9 




(A2» |s-t;i m) )_ 


- < 


II -R llg 






1(1+1) 

R 2 




i? 2 




.R 2 



< 



R 



cR 



ll^lloo|K m) L+Plloo||^i m) | 



< 



II^IL + cVaF 
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For the second set of inequalities: 



\(yim\v^) 



' X^|A^)>|<-^||^lU|^ m) || ff < 

An 



, (m) 
An 



c\\dy lm \\ d 



< 



A, t 



"0 



(m) 



For the last set of inequalities of lemma 4, once we have the results of the 
last proposition we can follow the way of [H, or that presented in |I. ■ 



3 Numerical application 

In general, the spectrum of the truncated matrices does not converge at the 
exact spectrum. Without additional results, one knows that only the lowest 
eigenvalue of the truncated matrices converges at the exact value. About 
these facts, one can consult ||. The results of the last section have another 
important consequence: in the proposed basis, the spectrum of the truncated 
matrices converges at the exact spectrum. Moreover, because the matrix 
of the Laplace-Beltrami operator in the y basis is " cvasidiagonal" in the 
sense that all nondiagonal elements are bounded by WhW^ and the diagonal 
elements increase approximatively as I (I + 1) /R 2 , it is to be expected that 
the spectrum of these truncated matrices to be very stable. That means, 
that even for low dimensions these matrices give us a good approximation of 
the exact spectrum. Let us choose the following particular surfaces for our 
numerical application: r (9) = 1 + 1.2 cos (6>) + 3 cos (#) 2 , presented in figure 
1. The eigenvalues for different truncated matrices and m = are presented 
in figures 2-6. 

Now, let us choose an orthonormal basis for which the affirmation of 
Lemma 4 is not true. If d\x g (9, <p) — a (9, (f>) sin 9 ■ d9d(j) is the measure 
induced by the metric g in the coordinates {9, 0}, then: 

Proposition 7 The set y of functions: 

yim (0,0) = -^==> m = 0,1,...,/= \m\ , \m\ + 1,..., 



11 



is an orthonormal basis in L 2 (M, /i g ). 

The proof of this proposition is analogous with that of Proposition 1. The 
eigenvalues of different truncated matrices, calculated in this basis and for the 
case m = 0, are presented in figures 7-11. The numerical application shows 
that in this case the spectrum of the truncated matrices are very unstable. 

This instability can be considered as an indicator of the fact that for the y 
basis the affirmation of our theorem is not true. 

4 Conclusion 

This paper has shown how to construct an orthonormal basis in the space 
of square integrable functions defined on a C°° surfaces with axial symme- 
try, basis which is appropriate for the problems which involve the Laplace- 
Beltrami operator. The procedure is standard, in the sense that it can be 
applied following the same steps for any C°° surface with axial symmetry. 
The stability of the truncated matrices spectrum was theoretically antici- 
pated and numerically verified. By practical point of view, this allow us to 
use truncated matrices with small number of rows and columns. 
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5 List of figures 

Fig.l The particular surface chosen for our numerical application: r (9) = 
1 + 1.2 cos (0) +3cos(#) 2 . 

Fig. 2 The eigenvalues of the 15x15 truncated matrix of Laplace-Beltrami 
operator in y orthonormal basis. 

Fig. 3 The eigenvalues of the 20x20 truncated matrix of Laplace-Beltrami 
operator in y orthonormal basis. 

Fig. 4 The eigenvalues of the 25x25 truncated matrix of Laplace-Beltrami 
operator in y orthonormal basis. 

Fig. 5 The eigenvalues of the 30x30 truncated matrix of Laplace-Beltrami 
operator in y orthonormal basis. 
Fig. 6 The superposition of figures 2-5. 

Fig. 7 The eigenvalues of the 15x15 truncated matrix of Laplace-Beltrami 
operator in y orthonormal basis. 

Fig. 8 The eigenvalues of the 20x20 truncated matrix of Laplace-Beltrami 
operator in y orthonormal basis. 

Fig. 9 The eigenvalues of the 25x25 truncated matrix of Laplace-Beltrami 
operator in y orthonormal basis. 

Fig. 10 The eigenvalues of the 30x30 truncated matrix of Laplace-Beltrami 

operator in y orthonormal basis. 

Fig. 11 The superposition of figures 7-10. 

Note: the horizontal coordinate in figures 2-11 is just an ordering index which 
puts the eigenvalues of the truncated matrices in an increasing order. 
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